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Abstract 

The shear free condition is studied for dissipative relativistic self-gravitating flu- 
ids in the quasi-static approximation. It is shown that, in the Newtonian limit, such 
condition implies the linear homology law for the velocity of a fluid element, only 
if homology conditions are further imposed on the temperature and the emission 
rate. It is also shown that the shear-free plus the homogeneous expansion rate condi- 
tions are equivalent (in the Newtonian limit) to the homology conditions. Deviations 
from homology and their prospective applications to some astrophysical scenarios 
are discussed, and a model is worked out. 
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1 Introduction 



As it is well known the shear plays an important role in general relativistic and cosmo- 
logical models (see Collins & Wainwright 1983 and Glass 1979 and references therein). 

In the case of slowly evolving (quasi-static) non-dissipative systems, it can be shown 
that in the Newtonian limit, the shear-free condition leads to the homologous contraction 
(or expansion) law for the velocity (Herrera & Santos 1995). However this is not neces- 
sarily the case in the presence of a heat flow vector (Herrera & Di Prisco 1996) and/or 
free streaming radiation (see below). This fact, and the great relevance of homology 
conditions in astrophysics (Schwarzschild 1958; Kippenhahn & Weigert 1990; Hansen & 
Kawaler 1994) provide the main motivation for this work. 

It is our purpose here to explore deeper the link between these two conditions and 
present some astrophysical scenarios where departures from homologous evolution (keep- 
ing the shear-free condition) migth drastically change the whole picture of the system. 

Accordingly, we shall consider dissipative systems. Indeed, dissipation due to the 
emission of massless particles (photons and/or neutrinos) is a characteristic process in 
the evolution of massive stars. In fact, it seems that the only plausible mechanism to 
carry away the bulk of the binding energy of the collapsing star, leading to a neutron star 
or black hole is neutrino emission (Kazanas & Schramm 1979). 

In the diffusion approximation, it is assumed that the energy flux of radiation (as that 
of thermal conduction) is proportional to the gradient of temperature. This assumption is 
in general very sensible, since the mean free path of particles responsible for the propaga- 
tion of energy in stellar interiors is usually very small as compared with the typical length 
of the object. Thus, for a main sequence star as the sun, the mean free path of photons 
at the centre, is of the order of 2 cm. Also, the mean free path of trapped neutrinos in 
compact cores of densities about 10 12 g. cm.~ 3 becomes smaller than the size of the stellar 
core (Arnett 1977; Kazanas 1978). 

Furthermore, the observational data collected from supernovae 1987A indicates that 
the regime of radiation transport prevailing during the emission process, is closer to the 
diffusion approximation than to the streaming out limit (Lattimer 1988). 

However in many other circumstances, the mean free path of particles transporting 
energy may be large enough as to justify the free streaming approximation. Therefore 
we will include simultaneously both limiting cases of radiative transport (diffusion and 
streaming out), allowing for describing a wide range situations. 

As mentioned before homologous evolution, an assumption widely used in astrophysics 
(Schwarzschild 1958; Kippenhahn & Weigert 1990; Hansen & Kawaler 1994), is known to 
be equivalent, in the non-dissipative case, to the shear-free condition in the Newtonian 
limit (Herrera & Santos 1995). As we shall see here, the presence of dissipative terms 
requires further the assumption of homology conditions on temperature and emission rate, 
in order to keep the homologous linear dependence of the velocity. It will also be shown 
that imposing the rate of expansion to be independent of the radial coordinate (together 
with shear-free condition) amounts to the full set of homology conditions. 

Although deviations from the homologous evolution are shown to introduce extremely 
small modifications in the expression for the velocity, these terms might be relevant in 
some very specific situations which we will discuss later. 
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It is also worth mentioning that although the most common method of solving Ein- 
stein's equations is to use comoving coordinates (e.g. May & White 1966; Wilson 1971; 
Burrows & Lattimer 1986; Adams, Cary & Cohen 1989; Bonnor & Knutsen 1993) we 
shall use noncomoving coordinates, which implies that the velocity of any fluid element 
(defined with respect to a conveniently chosen set of observers) has to be considered as a 
relevant physical variable (Bonnor & Knutsen 1993). 

The plan of the paper is as follows. In Section 2 we define the conventions and give 
the field equations and expressions for the kinematical and physical variables we shall use, 
in noncomoving coordinates. In Section 3 we give the general expression for the velocity 
and evaluate the dissipative terms. A very simple model is presented in Section 4. Finally 
a discussion of results is presented in Section 5. 



2 Relevant Equations and Conventions 
2.1 The field equations 

We consider spherically symmetric distributions of collapsing fluid, which for sake of 
completeness we assume to be locally anisotropic, undergoing dissipation in the form of 
heat flow and/or free streaming radiation, bounded by a spherical surface E. 
The line element is given in Schwarzschild-like coordinates by 

ds 2 = e v dt 2 - e x dr 2 - r 2 (d6 2 + sin 2 6d<p 2 ) , (1) 

where v(t,r) and A(t, r) are functions of their arguments. We number the coordinates: 
x° = t; x 1 = r; x 2 = 8; x 3 = 0. The metric has to satisfy Einstein field equations 

g; = -sttt;. (2) 

In order to give physical significance to the T£ components we apply the Bondi approach 
(Bondi 1964). Thus, following Bondi, let us introduce purely locally Minkowski coordi- 
nates (r, x, y, z) 

dr = e u ^ 2 dt ; dx = e x ^ 2 dr ; dy = rdO ; dz = r sin 9d(f). 

Then, denoting the Minkowski components of the energy tensor by a bar, we have 

f ° = T ° ; fl = Tl ; f| = T 2 ; f 3 3 = T 3 3 ; f 01 = e~ {u+x)/2 T 01 . 

Next, we suppose that when viewed by an observer moving relative to these coordinates 
with proper velocity uj in the radial direction, the physical content of space consists of 
an anisotropic fluid of energy density p, radial pressure P r , tangential pressure P±, radial 
heat flux q and unpolarized radiation of energy density e traveling in the radial direction. 
Thus, when viewed by this moving observer the covariant tensor in Minkowski coordinates 
is 



P + e 


-q-e 





-q-e 


P r + e 








P± 











Pa 
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Then a Lorentz transformation readily shows that 



_ - _ p + P r u 2 2Quoe x / 2 

T - T - 1_ U 2 + (1-^2)1/2 + 6 ' ^ 

(1-^)1/2 6 ' ^ 

T* = T* = T* = T* = -P ± , (5) 



T = e (^+A)/2y 



01 



(p + P r )ue^ +X ^ 2 
1-u 2 



2 \ _ ^("+>)/2 



(1 +lu 2 ) -ee 



(6) 



with 



Q 



qe 



-A/2 



and 



(l-cu 2 ) 1 ^ 
1 + u 



(7) 



e = e- 



1-u/ 



Note that the coordinate velocity in the (t, r, 9, 0) system, dr/dt, is related to uj by 



uj = -e^ 2 . 
dt 



Feeding back (jHHEl) into (J2J), we get the field equations in the form 



p + P r uj 2 2Que x ' 2 _ 1 

l-u 2 + (i-u 2 y/ 2 + e ~~8^ 

P r + puj 2 2Qooe x / 2 _ 1 

l-u 2 + (i -u 2 y/ 2 + e ~ 



1 



1 i/ N 



(9) 

(10) 
(11) 



Pi 



1 f e~ l 
.11 



2v" + v' 2 - AV + 2 



2A + A(A-z>) 
z/-A'V 



(12) 



(P + fiW A)/2 + Q^t^jl + u 2 ) + ee^ 2 



A 



(13) 



1 - UJ 2 (1 - Co' 2 ) 1 /2 • 87rr 

where the dots and primes stand for partial derivatives with respect to t nd r respectively. 
Outside of the fluid distribution, the spacetime is that of Vaidya, given by 
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ds 2 



1 - 



2M(u) 
R 



du 2 + 2dudR - R 2 (d6 2 + sin 2 8d<f) 2 ) , 



(14) 



where u is a coordinate related to the retarded time, such that u = constant is (asymp- 
totically) a null cone open to the future and R is a null coordinate (gRR = 0). It should be 
remarked, however, that strictly speaking, the radiation can be considered in radial free 
streaming only at radial infinity. The two coordinate systems (t, r, 8, 0) and (u, R, 8, 4>) 
are related at the boundary surface and outside it by 



u = t — r — 2M In ( ^— - 1 
\2M 



(15) 



R = r. (16) 

In order to match smoothly the two metrics above on the boundary surface r = r 2 (t), we 
first require the continuity of the first fundamental form across that surface. Which in 
our notation implies (Herrera et al. 2002) 



= 1 



= 1 



2M 
2M 



Ry 



;i7) 



(18) 



Where, from now on, subscript S indicates that the quantity is evaluated at the bound- 
ary surface £ and R = R-e(u) is the equation of the boundary surface in (u,R,0,4>) 
coordinates. And 



[Prk= Qe X / 2 (l-LU 2 ) 



1/2 



(19) 



expressing the discontinuity of the radial pressure in the presence of heat flow, which is a 
well known result (Santos 1985). 

Next, it will be useful to calculate the radial component of the conservation law 



T M = 0. 



After tedious but simple calculations we get 

16"7T 



(7? - T|) + W (Tl - T °) + ^ (A + ^ - ^ 



(20) 



(21) 



which in the static case becomes 



v' , ^ x 2 (P ± - P r ) 



(22) 



representing the generalization of the Tolman-Oppenheimer-Volkof equation for anisotropic 
fluids (Bowers & Liang 1974). 
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2.2 The kinematical variables 

The components of the shear tensor are defined by 



&flV ~\- Uy^ U^CLy UyCL^ Q P ^y , 



(23) 
(24) 



where 

denote, P^ u the projector onto the three space orthogonal to w M , 6 the expansion and a M 
the acceleration. 

A simple calculation gives 
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,-f/2 



2(l-cu 2 ) 1/2 V 1-u 



A + 



2uouj 



e- A /2 / 2c; V 4uA 

+ — I-T/o ( ^ + 2u/ + ; + — , 



2(l-cu 2 ) 1/2 



1 — a; 2 r / 



(25) 



3 (1 - cu 2 ) 3/2 

+e V2 ( Wl/ + 



2a; ci; 
2uj' 2u 



1-u 2 



e -\ 2 (l-u 2 ) 

(722 = " fll, 



0"33 



e -A r 2 (1 _ ^2) 



sin 2 9a u , 



(26) 

(27) 
(28) 



(29) 



001 



-w e ("- A )/Vii, 



and for the shear scalar, cr = (cr flv a fJ ' v ) 1 / 2 , 



a 



V3 



e e - A / 2 

■3 ~(1 -cu 2 )V2 



2.3 The Weyl tensor 



(30) 



(31) 



The model to be presented in Section 4 is obtained from the assumption of conformal 
flatness. Furthermore, since the publication of Penrose's (1979) work, there has been an 
increasing interest in studying the possible role of Weyl tensor (or some function of it) in 
the evolution of self-gravitating systems (Wainwright 1984; Goode & Wainwright 1985; 
Bonnor 1985, 1987; Goode, Coley & Wainwright 1992; Pelavas k Lake 2000 ; Herrera et al. 



7 



2001). This interest is reinforced by the fact that for spherically symmetric distribution 
of fluid, the Weyl tensor may be defined exclusively in terms of the density contrast and 
the local anisotropy of the pressure (see below), which in turn are known to affect the fate 
of gravitational collapse (Mena & Tavakol 1999; Eardley & Smarr 1979; Christodoulou 
1984; Newman 1986; Waugh & Lake 1988; Dwivedi & Joshi 1992; Singh & Joshi 1996; 
Herrera & Santos 1997; Bondi 1993; Barreto 1993; Coley & Tupper 1994; Martinez, Pavon 
& Nunez 1994; Singh, Singh & Helmi 1995; Das, Tariq & Biech 1995; Maartens, Maharaj 
& Tupper 1995; Das et al. 1997; Corchero 1998a,b; Bondi 1999; Hernandez, Nunez & 
Percoco 1999; Harko & Mak 2000; Das & Kloster 2000; Joshi, Dadhich & Maartens 2002; 
Krisch & Glass 2002; Corchero 2002; Harko & Mak 2002; Mak & Harko 2002). 

Therefore it is worthwhile to include here some expressions for the Weyl tensor. Thus, 
using Maple V, it is found that all non-vanishing components of the Weyl tensor are 
proportional to 

W = T -C! 32 = W {8) + r -^f (A + £ - *fj , (32) 

where 

e A 1 v'X v a v" X' u'\ 

^~^ + ~"T~y~27 + 27j' (33) 

corresponds to the contribution in the static (and quasi-static) case. Also, from the field 
equations and the definition of the Weyl tensor it can be easily shown that (see Herrera 
et al. 1998 for details) 

W = ~ [ r 3 (T )' dr + ^r 3 (t| - I*) . (34) 
2.4 The slowly evolving approximation 

In this work we shall consider exclusively slowly evolving systems. That means that our 
sphere changes slowly on a time scale that is very long compared to the typical time in 
which it reacts on a slight perturbation of hydrostatic equilibrium, this typical time is 
called hydrostatic time scale. Thus our system is always in hydrostatic equilibrium (very 
close to) and its evolution may be regarded as a sequence of static models linked by (fT5j) . 

This assumption is very sensible because the hydrostatic time scale is very small for 
almost any phase of the life of a star. It is of the order of 27 minutes for the sun, 4.5 
seconds for a white dwarf and 10~ 4 seconds for a neutron star of one solar mass and lOKm 
radius (Schwarzschild 1958; Kippenhahn & Weigert 1990; Hansen & Kawaler 1994). 

Let us now express this assumption through conditions for u and metric functions. 

First of all, slow contraction (or expansion) means that the radial velocity u measured 
by the Minkowski observer, as well as time derivatives are so small that their products as 
well as second time derivatives can be neglected. Thus we shall assume 

v « A w Az> w A 2 w z> 2 w u 2 w Co ta 0. (35) 



w M = ^ ( 
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Then, it follows from (|13|) that Q and e are, at most, of order 0(u). Henceforth, with 
this approximation, (J21j) becomes 

(P r + e)' + (p + P r + 2e) V - - 2 P± ~ Pr ~ e = (36) 

2 r 

which is the equation of hydrostatic equilibrium for an anisotropic fluid radiating a null 
fluid of energy density e. 

Thus, as mentioned before, the system, although evolving, is in hydrostatic equilibrium 
(up to order 0(u)), this allows for a very simple extension of any static solution to the 
slowly evolving case. 



3 Shear— free and homology conditions 

As mentioned before the only relevant component of the shear tensor is an given by 
equation (|26[). Evaluating this last equation in the slowly evolving approximation, we 
obtain 

-»/*X + e- x ' 2 (uv ' + 2u/-^)]. (37) 



2 A 



Next, using (|T3*jl and 

P r + P=^—(i/ + \ , )-2e, (3* 
easily obtained from (fTUj) and (fTT|) . one gets 

on = =-e = — , 2 e 

r z r z sin # 

= _l e A / 2 L' - — - - - 4vrrQe 3A / 2 - 47rree A ^ . (39) 
3 \ 2 r / 

We can solve for cu, to obtain 

w = e^ A_As ^ 2 - 4vrre A / 2 P ( Qe A + ee A / 2 - — e - A ^) dr. (40) 
r s Jr V 167r r / 

From the above equation we find that in the non-dissipative, shear-free case we obtain 



uj = ^—e^-^ 2 . (41) 

While in the Newtonian limit we have M(u) ~ A ~ v ~ and we recover the well known 
linear expression, typical of the homologous evolution (Schwarzschild 1958; Kippenhahn 
& Weigert 1990; Hansen & Kawaler 1994), 

r 

^Newt = W S ■ (42) 

Also, from ()31|) evaluated in the slowly evolving approximation, it follows that in the 
shear-free motion 

= ^e" A / 2 , (43) 
r 

9 



which of course is valid also in the dissipative case. Using (J41|) . we can write 



rs 



(44) 



implying that even in the general (relativistic) case, the expansion rate is homogeneous 
(independent of r) for the slow, and dissipativeless shear-free motion. 

Let us now consider the dissipative shear-free case. 

From the relativistic Maxwell- Fourier law, we have 



or 



Q 



Q = — Ke 



-A / rpl 



(45) 



(46) 



where T is the temperature and n denotes the coefficient of conduction. It should be 
reminded that in the quasi-static approximation, the system is assumed to be relaxed at 
all times (the relaxation time is zero) and accordingly, any hyperbolic transport equation 
reduces to (J43j) . 

Then feeding back (|47)|) into (|4T)|) and using (fT%|) together with the shear-free condition, 
we obtain 



to 



we— 1 



2M(u) 
rs 



1/2 



+ 47T/t (T s - T) r 



''E 



+2iiKr J TiAlr— Inr I ee A ^ 2 c?r 

which in the Newtonian limit yields 

r f 

u Newt = ws h 47r k (Ts -T)r — Anr / 

r s Jr 



D \/2 



''E 



edr. 



Also, it follows that the expansion (47) with (51) can be written 



e = e s + 3 



Attk (Ts - T) + 2ttk / Tv'dr -Air ee 

J r Jr 



(47) 



(48) 



(49) 



Thus, unlike the non-dissipative case (see also Herrera Sz Di Prisco 1996), the shear- free 
collapse in the Newtonian limit does not yield the linear law of homologous contraction 
(Schwarzschild 1958; Kippenhahn & Weigert 1990; Hansen & Kawaler 1994), unless we 
impose further homology conditions on T and e, i.e. unless we assume that for any given 
fluid element, all along the evolution 



TV 



constant, 



— = constant. 
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From (47) we observe that the sign of u for any value of r, is not necessarily the same as 
that of (as is the case in the non-dissipative evolution). In particular, for sufficiently 
large (negative) gradient of temperature and/or sufficiently large (positive) e term, we 
may have > and uj < 0. The same conclusion of course applies to 0. 

In other words the system may be evolving in such a way that inner shells collapse, 
whereas outer ones expand. 

This effect, which we have called "thermal peeling" (Herrera & Di Prisco 1996), is 
also present in the relativistic regime, provided the third term in the right side of ()47|) 
is not too large. It represents the analog of the "cracking", however whereas the later 
takes place, under some conditions, when the system abandons the state of equilibrium or 
quasi-equilibirum (Herrera 1992), the former occurs while the systems is evolving quasi- 
statically. 

However, observe that expressing variables in c.g.s. units, we have that, 

kT~1(T 59 \KWT\cm- 1 , 

where [«] and [T] denote the numerical values of these quantities as measured in 
er g s _1 arr x K^ 1 and K respectively. Therefore extremely high conductivities and/or 
AT are required for thermal peeling to be observed in the Newtonian regime. Also, we 
have 

e ~ 1(T 59 [e]cm- 2 , 

where [e] denotes the numerical values of this quantity as measured in erg s~ l cmT 2 . 

Before closing this section, it is worth mentioning that, in general, such high thermal 
conductivities are associated to highly compact, degenerate objects where the Newtonian 
limit is not reliable. 

Also, it should be noticed that in (}4~5j) it has been assumed that terms of order 0(M/r^) 
and higher are negligible with respect to k (T^ — T). This of course is not always true, 
as commented above, in which case (J48j) is not valid. Finally, it is worth noticing that 
demanding to be homogeneous, we are lead to the homologous contraction, implying 
thereby that (in the Newtonian limit), the shear-free and homogeneous expansion rate 
conditions are equivalent to the whole set of homologous conditions. 



4 A model 

In order to illustrate the point raised in Section 3, let us present a very simple model 
based on the assumption of conformal flatness and shear-free condition. Also, since local 
anisotropy does not enter explicitly in ()47|) we shall assume P r = P±. 
Thus assuming W — 0, it follows from 



ft = -• (50) 
r 



Next, taking for simplicity Q = (pure free streaming dissipation) and 



r 



(3r ( 1 - - ) , (51) 
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with (3 = j3(t) one obtains 

p = 3/3r(l-^)+ 7 (t), (52) 

where 7 = p(0,t). Observe that with this choice of e (if we assume Q — 0, i.e the 
dissipation takes place at the free streaming approximation exclusively), the evolution 
proceeds adiabatically (the total mass is constant) even though e 7^ within the sphere. 
Next, from the definition of the mass function (Herrera & Santos 1995) 

m(r, t) = 4vr f r 2 T °dr = Ul - e~ A ) (53) 
Jo 2 

and junction conditions, it follows 

m = — (3^ + 7--^— J , (54) 

and 

2 / 3M \ 

" = ^(s^ -'J- (55) 

with M = mj. As expected, if we put /? = we recover the well known interior 
Scwarzshild solution (evolving quasi-statically). The remaining of the metric and physical 
variables may now be easily obtained from the field equations. Feeding back (|51|) into 
(J4*U)) one sees that playing with (3 it is possible (at least in principle) to obtain u < for 
some values of r, even though u>^ is assumed possitive (peeling). 



5 Conclusions 

We have seen how dissipative terms affects the radial dependence of to and the expansion 
rate, in the shear-free case, if we relax the homology conditions on dissipative variables. 
We have also seen that the dissipative terms may lead to a "peeling". However these 
contributions appear to be extremely small and therefore it is pertinent to ask if there 
exist astrophysical scenarios where dissipative contributions might have some effect on u, 
and in particular if they could produce a "peeling" . 

Assuming the highest values for luminosity at the last stages of stellar evolution, of 
the order of 10 5 times the sun luminosity, produced at a shell of radius of 1/10 of solar 
radius, we only get 

e ~ 1(T 36 cm' 2 . 

A more promising case is provided by the Kelvin-Helmholtz phase of the birth of a 
neutron star (Burrows & Lattimer 1986). Indeed in this phase, during tens of seconds, 
some 10 53 ergs are radiated away. If this energy is transported via diffusion to the surface, 
then assuming (Flowers & Itoh 1979,1981) 

k w 10 23 [p/10 u g cm~ 3 ][10 8 K/T]erg s' 1 cm' 1 K" 1 , (56) 

we see that the corresponding contribution to (|47|) is still too small. However if we assume 
that part of 10 53 ergs are propagated in the free streaming regime, then the last term in 
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(|47|) for sufficiently small r (as compared to r^) is of the order of 10 52 /r£. Therefore for 
positive surface velocities of the order of 30m/ s there may be a peeling (u < for r < r^). 

Finally, let us mention that in a pre-super novae event, values of the order of 10 13 and 
10 37 have been estimated for [T] and [k] respectively (Martinez 1996). With these values, 
it is clear that peeling is also possible, in particular for sufficiently large values of r^. 
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